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Abstract 

Thanks to the unnaturally small value of the QCD vacuum angle 0 < 10“^®, time- 
reversal violation (2^ offers a window into physics beyond the Standard Model (SM) 
of particle physics. We review the effective-field-theory framework that establishes a 
clean connection between^mechanisms, which can be represented by higher-dimensional 
operators involving SM fields and symmetries, and hadronic interactions, which allow for 
controlled calculations of low-energy observables involving strong interactions. The chiral 
properties ofp mechanisms leads to a pattern that should be identifiable in measurements 
of the electric dipole moments of the nucleon and light nuclei. 
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1 Introduction 

Violation of time reversal (T) is a fundamental asymmetry between past and future, the 
microscopic dynamics not being invariant under change in the sign of time. In a Lorentz- 
invariant quantum held theory, where the product CPT is conserved, T violation is 
equivalent to violation of the product of charge conjugation (C) —the exchange between 
particle and antiparticle— and parity (P) —the change of sign in spatial coordinates. CP 
violation (jCP^ is one of the ingredients needed to explain why the visible universe seems 
to be made predominantly of matter, without a signihcant fraction of antimatter [^. 

The Standard Model (SM) contains a source of jVP and^: the phase of the CKM 
matrix [^, which appears in observables through a combination J ~ 3-10“^ of matrix elements 
[^. While this mechanism explains the violation observed in K and B decays (^, it gives only 
very small contributions to quantities that do not involve flavor change between initial and hnal 
states. In particular, it is not sufficiently large to account for the observed matter-antimatter 
asymmetry [^. 

The ideal observables to probe new^ interactions are flavor-conserving observables such as 
^electromagnetic form factors (EEs), which can be split into electric, magnetic, and toroidal, 
depending on whether they interact with long-range electric, long-range magnetic, or short- 
range electromagnetic helds. Their P and T transformation properties are summarized in 
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Multipolarity 

Electric 

Magnetic 

Toroidal 

0 (monopole) 

PT 

— 

— 

1 (dipole) 


PT 

/T 

2 (quadrupole) 

PT 


P/ 


Table 1: Parity and time-reversal properties of electromagnetic form factors according to 
multipolarity. S {0) denotes that the symmetry S is preserved (violated). The pattern 
repeats as multipolarity increases. A particle of spin s has multipoles up to 2s. 


TableOf particular interest are permanent electric dipole moments (EDMs), which require 
both ^and^. Even with current technology they effectively probe very small distances. For 
example, the existing bound on the neutron EDM, \dn\ < 2.9 ■ 10“^^e fm [^, means that a 
charge imbalance, if any, effectively takes place at a distance 13 (or more) orders of magnitude 
smaller than the size of the neutron. Other ^ multipoles, like magnetic (MQM) and toroidal 
(TQM) quadrupole moments, are less accessible experimentally. 

Searches are in progress around the world for EDMs of the neutron and of (neutral) atoms 
and molecules, which are sensitive to the EDMs of the electron and nuclei, and their ^ inter¬ 
actions. A new generation of experiments 10 promises to improve neutron EDM sensitivity 


by one or two orders of magnitude, which is remarkable but still above the expected CKM 
“background” at ~ 10“^®e fm (see Ref. for an assessment and references to original papers, 
for a recent discussion). Even more exciting is the groundbreaking proposal (see 


and Ref. 

12 

Refs. 13 

14 


13,14 for summaries and references) that the EDMs of charged particles be investigated 


in specihcally designed storage rings, and not just as byproducts of other experiments as for 
the muon [^. We might see the deuteron EDM (d^) probed at the level of ~ 10“^®e fm [Id] , 
and similarly for the EDMs of the proton (dp) and of the nucleus of helium-3, helion (d/i). 

The discovery of an EDM above the CKM background would be a signal of long-sought 
new physics, but would, by itself, leave us in the dark about its origins. From a theoretical 
perspective, it is important to investigate the set of EDMs that would allow us to identify the 
dominant source(s) of T- Such an identihcation is likely to give clues about physics beyond 
the SM (BSM) and the scale of this new physics, which we denote Mf. The aim of this 
review is to show that the framework of effective held theories (EFTs), coupled to recent and 
not-too-distant-future progress in strong-interaction physics, will allow us to carry out this 
identihcation for the T operators that involve quarks and gluons. 

There are good experimental and theoretical reasons to believe that the SM is an EFT for 
processes involving momenta Q ~ Mew ~ 100 GeV. In addition to the possible existence of 
new light or stable heavy particles, new physics can be represented at the electroweak (EW) 
scale Mew by operators of canonical dimensions d > 4, which involve known particles and are 
constrained by the SM symmetries, namely Lorentz invariance and gauged SU{3)c x SU{2)i x 
U{1)y. One expects these operators to have strengths 0{M^ '^), making the lowest-dimension 
operators most signihcant. 

The SM has other d = 4 ^ operators, which involve the non-Abelian gauge bosons and, 
despite being total derivatives, can contribute to observables due to topological ehects 
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The operator involving gluons could give rise to large EDMs of hadrons and nuclei, 
but the neutron EDM bound already signihcantly constrains its dimensionless strength, the 
QCD vacuum angle This unnatural value is the famous “strong CP problem”. 

The most promising solution is offered by the Peccei-Quinn (PQ) mechanism 17 , where 


an additional approximate symmetry, 17(1)pq, is spontaneously broken generating a small 
6 dynamically. The corresponding pseudo-Goldstone boson, the axion 18,19 , is a viable 


dark-matter candidate. Whatever the mechanism may be, the smallness of 9 leaves room for 
higher-dimensional operators. 

The sole d = 5 interaction among known particles 


20 gives rise to neutrino masses and 


lepton-number violating processes, which are searched for with neutrinoless double-beta decay 
(see Ref. 21 for a recent review). One can expect violation from phases in the correspond¬ 


ing PMNS mixing matrix [^, but it is unclear if the observed baryon-antibaryon asymmetry 
can be generated through leptogenesis 


23 , a mechanism based on the simplest ultraviolet 


(UV) completion of the d = 5 interaction. 
The many d = 6 interactions 


24 have been conveniently cataloged in Ref. 25 . Among 


those responsible for hadronic and nuclear EDMs are most sensitive to quark EDMs 

(qEDMs) and chromo-EDMs (qCEDMs), the gluon chromo-EDM (gCEDM) [^, and certain 
four-quark interactions 28,^. Nowadays the techniques exist—renormalization-group (RG) 
running down to the QGD scale Mqcd ~ 1 GeV, lattice QGD (LQGD) for the calculation of 
low-energy constants (LEGs), nuclear EFTs to describe the dynamics at momenta Q < Mqcd 
in terms of the LEGs— to connect these operators at the EW scale to light-nuclear EDMs. 

A frequent misconception is that nuclear-physics errors would obfuscate any of the minute^ 
effects we are interested in. This presumption is obviously not correct for quantities like EDMs, 
which vanish when ^ parameters vanish; in this case errors affect only the proportionality 
factor, and as we are going to see they should not affect many of the conclusions. Light nuclei 
are special in this regard. A key ingredient is the approximate chiral symmetry of QGD, and 
its breaking. Not only does chiral symmetry provide the basis for a systematic expansion of 
hadronic and nuclear observables in powers of Q/Mqcd 30 , but it also acts as a “hlter” to 
separate the effects of various ^ sources. While all these sources generate ^ observables, they 
break chiral symmetry in different ways and as a consequence produce different patterns in 
the relative magnitudes for these observables. Ghiral symmetry was already an important 
aspect of the classic studies of the neutron EDM from 6 31,32, and now it has been extended 
to other sources and more nucleons 33 -36 . A single measurement (say, dn) can always be 


attributed to any one source (say, a 0 of just the right, minute size), but, as we review in 
Sections and 1^ combined measurements of dn, dp, dd and dh provide increasingly detailed 
information on the ^ sources. The further measurement of the triton EDM (dt) would allow 
as good a separation of underlying mechanisms as possible in the strong-interacting sector 
at low energies, under the assumption that lower-dimension operators are most important. 
(Measurements of the deuteron MQM and TQM would also be valuable, but seem impossible 
for the foreseeable future.) 

All the main techniques needed for this analysis have experienced signihcant progress 
recently. RG running down to the QGD scale has long been used as a tool to investigate the 
low-energy consequences of specihc BSM models, and a model-independent summary of these 
results has appeared recently 37 . Gomputational advances are bringing LQGD to the forefront 
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of hadronic and nuclear physics, and the time is approaching when EFTs will be able to use 
LQCD data, rather than experiment, as input in the calculation of nuclear properties 38 


Still, the study of 'X matrix elements is in its infancy, particularly for d > 4 operators, and 
this constitutes the biggest gap in connecting nucleon and light-nuclear EDMs to BSM ^ 
interactions. Nuclear EFTs make it possible to approach hadronic and nuclear physics 
incorporating SM symmetries, and the ^ hadronic interactions involve, in lowest order, six 
LECs. Nuclear potentials inspired by Chiral EFT 39,40 , where pions and chiral symmetry 
play a signihcant role, are now the favorite starting point for “a& initio" nuclear-structure 
methods. Although fully consistent calculations are not yet possible, for light nuclei several 
tests suggest that errors are relatively small and can be better quantihed once some subtleties 
in the RG of pion exchange are clarihed 41,42. A signihcant future step would be to extend 


the framework summarized here to heavier nuclei, in order to enable a consistent analysis of 
atomic/molecular experiments as well. 

We limit ourselves here to a review of the techniques that address ^ in the nucleon and 
light nuclei from a model-independent perspective. The relevant operators at the quark/gluon 
level, including the RG to Mqcd, are introduced in Section and translated into hadronic 
interactions in Section The PQ mechanism is briehy summarized in Section Sections 
and 1^ review calculations of T electromagnetic observables for the nucleon and light nuclei, 
respectively, starting from the hadronic interactions. An outlook is reserved for Section 
A much more comprehensive review —which includes implications of specihc BSM models, 
calculations of QCD matrix elements with various assumptions, calculations of heavy-nuclear 
T quantities, a discussion of atomic/molecular EDMs, and many more references to earlier 
work— has appeared recently |^, and we refer the reader to it for a more detailed look at 
how the program presented here relates to other efforts. 


2 T the Quark-Gluon Level 

We summarize here the most important T interactions among quarks, gluons and photons. 
This background will allow the construction of T' hadronic interactions in Section 

2.1 QCD 0 Term and Higher-Dimensional Operators 

BSM physics can be described below its characteristic scale Mf > Mew in terms of SU (3)c x 
SU{2)l X U{1)y and Lorentz-symmetric operators. The kinetic terms involving quarks inter¬ 
acting with gluons (a = 1, ...,8) and weak bosons {i = 1,2,3) and with strengths 
( 7 s, g, and g', respectively, are 

= qLipqL + URipUR + dmpdR - i , (1) 

where q^ is a doublet of left-handed quarks; ur and dR are right-handed up- and down- 
type quarks; — igsG‘l^ta — igW^ti — ig'B^Y is the gauge covariant derivative with 

ta = Aa/2 (Aa are the Gell-Mann matrices), = rj/2,0,0 (r* are the Pauli matrices) and 
Y = 1/6,2/3,—1/3 for qi, ur and dR, respectively; and = dfj,Gl — d^G^ — gsf°'^'^G^^Gl, 
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and = d^B^ - d^B^ are, respectively, the SU{3)c, 
SU ( 2 ) £, and U{1)y held strengths with strnctnre constants for SU ( 3 )c and e*-^^ for SU ( 2 ) l. 

For simplicity we omit indices that rnn throngh the three generations, which are snmmed over. 
The d = 4 terms related to ^ are the Ynkawa conplings of the qnarks and the topological 


9 term 16 


cf = - UR + dn) + H.C. - , 


( 2 ) 


where (f is the Higgs donblet and (Here = +1.) The Ynkawa 

conplings form 3x3 complex matrices in flavor space. Relative phases lead to in the 
CKM matrix and, as described below, there is an interplay between the overall phase and 
the vacnnm angle 0 < d < 27r. A term analogons to 9 for weak gange bosons gives negligible 
contribntions at low energies. 


Becanse of the smallness of d = 4 we follow Refs. 35,36 and consider the d = 6 


terms 25 


L 


( 6 ) 

r 


= -2 




{qLY'^pun + qLY'^pdn) + H.c. + 


[ 9 sT^taG% + + g'TlB,, )puR + H.c. 


( QstXGl, + gVi^nWy + )pdR + H.c. 


+ gl p^iD.p + H.c. 

+gl [Si qiuR q^dR + Eg qitaUR gftadij] + H.c., 
where v ~ 246 GeV is the Higgs vacnnm expectation valne (vev). Here: 


( 3 ) 


The 3 X 3 complex matrices and the angle 9' correct, after EW symmetry breaking, 
Ynkawa conplings and 9 term. Terms linear in the Higgs held are of more phenomeno¬ 
logical interest, as they either modify the Ynkawa conplings of the Higgs to qnarks 
introdncing, in general, havor-changing effects, or affect the glnon-fnsion prodnction 


mechanism at the LHC 44 -47 


The coefficients and T are 3x3 matrices in havor space. The nondiagonal entries 


B,W 


contribnte to havor-changing cnrrents and play an important role in havor physics 48 


[4^ . The diagonal components determine the qnark electric (qEDM) and chromoelectric 
(qCEDM) dipole moments, as well as the qnark weak EDM, which plays only a minor 
role at low energy. 


The parameter dw of the Weinberg three-glnon operator 27 can be thonght of as the 


glnon chromoelectric dipole moment (gCEDM). Similar terms involving weak gange 
bosons (for example, the boson’s EDM and MQM [^) can be relevant at colliders, 
bnt are very small at energies below Mew 
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• The Si s are complex four-index tensors in flavor space. If one considers only quarks of 
the hrst generation, there are two four-quark operators that respect SU{2)i (^ , 
which we refer to as PS4QOs due to their pseudoscalar nature. Additional four-quark 
operators can be constructed involving quarks of different generations. 


• The complex 3x3 matrix Si couples bosons to the right-handed quark current. 
It leads below Mew to additional four-quark operators , which couple left- and 
right-handed quarks and we call LR4QOs. It also causes in nuclear /3 decay when 
connect to the left-handed lepton current 

The d = 6 coefficients depend on the spectrum and parameters of the BSM model of 
choice, and can be determined by matching just below the scale M^. While a detailed study is 
beyond the scope of this review (see instead Ref. |^), from our brief discussion it is apparent 
that there is a rich interplay among the constraints on the couplings in Eq. ([^ that can be 
extracted from collider, flavor, and low-energy precision experiments. 

Below the scale Mew fhe breaking of SU{3)c x SU{2)i x U{1)y to SU{3)c x U{l)em is 
important and generates masses for fermions and weak gauge bosons. For studying low- 
energy observables like EDMs, the Lagrangian ([^ at the EW scale needs to be matched onto 
a theory with only light quarks, gluons and the photon = sin 9wW^ + cos where 

gsinOw = g' cos 6w = — e < 0. For nuclear physics applications we can limit ourselves to two 
light flavors, and at Mqcd the kinetic terms can be written as 

= iLiiPqL + inmn - (4) 
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where q = {u d)^, — igCJ^ta — ieA^Q is the SU{3)c x U{l)em covariant derivative 

in terms of the charge matrix Q = 1/6 -f tz/2, and F^i, = — d^A^ is the U{l)em held 


strength. Some of the issues associated with strangeness are discussed in Section |5.3 

from the QCD 6 term is intimately related to the quark masses. All the phases of 
the quark mass matrix can be eliminated through non-anomalous SU{2) vector and axial 
rotations, except for a common phase p, leaving 


4"^ = - {e^^qLMqn + e-^'^qnMqL) - 9 


_9^s_ 111/a^ f^a 

647^2 


(5) 


where M = m(l — ers) is the diagonal quark mass matrix with fa (the average light-quark 
mass) and e (the relative light-quark mass splitting) real parameters. The phase p and the 9 
angle are not independent. By performing an anomalous axial U{1)a rotation, all (AP can be 
rotated into the 9 term or into the complex mass term, and physical observables depend only 
on the combination 9 = 9 + 2p. For our discussion in Chiral EFT, it is convenient to eliminate 
the 9 term in favor of a complex quark mass. The additional SU{2) approximate symmetry 
of the Lagrangian can be exploited to align the vacuum in the presence of jAP to the usual 
QCD vacuum 31,32. From a low-energy point of view, vacuum alignment is equivalent to 


setting to zero the coupling of the neutral pion to the vacuum, at lowest order in the Chiral 


EFT expansion 33 . The Lagrangian becomes 


4^ = —mr(6*) qq + mr ^{9) q iers-l --—sind | q, 


1-e^ 


( 6 ) 
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where r{x) is an even fnnction of x, 


r(x) = 


/l + e^tan^I 


= 1 - n 






( 7 ) 


The 0(0^) terms become important in the PQ mechanism (Section]^. The last term in Eq. 
([^ is approximately linear in 9, and responsible for 

In the absence of flavor change, the d = 6 ^ operators that receive tree-level contribntions 
and are not snppressed by powers of Mew are 

^^6) ^ ^ ck (^dg + dsTs^ 

.9yV rabc fiuafi/^a /^b /^cp 
-r-^ysj t '^ap'^pp'^v 

_ 

+ ^ [ImSi {qqqi'f^q - qrq ■ qri'y^q) + ImSs {qtaqqilbtaq - qrtaq ■ qTi'ybtaq)] 

g2 ,, 

+ y (ImHi qnYq q'Tj'^plbq + ImHg qn^taq qrjlplbtaq) ■ (8) 

The tree-level matching of Eq. ([^ to Eq. ([^ is, for most operators, trivial. The qEDMs and 
qCEDMs are 


0,3 = -^Im [(T 


m (in.3 = 


+ rw)iiT2 


(r| - r^) 


iJ , md, 


0,3 


= -Im 
2 


r“±r 
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(9) 


where the indices refer to generations. In most models and (f)ii are proportional to 

the light-qnark Ynkawa conplings, thns canceling the light-qnark mass on the left-hand-side of 
Eq. ([^. The gCEDM gets tree-level contribntions only from itself, and similarly for PS4QOs. 
The imaginary part of the right-handed cnrrent in Eq. (|^ contribntes to the S'f/(2)L-breaking 
LR4QOs. Thns, at tree level. 


ImSi^s = , ImHg = 0, (10) 

where Vud — 0.97 is the np-down CKM element. ImSg is, however, generated by the QCD 
evolntion, as we discnss next. Note that all interactions in Eq. ([^ are also^; interactions 
are effectively of higher order and expected to prodnce even smaller effects. 


2.2 Sizes and Runnings of Couplings 

As the scale /x of interest decreases, one has to evolve the QF^ coefficients. The 9 term is not 
mnltiplicatively renormalized becanse 9 is periodic, althongh it can mix with the divergence 

The operators in Eq. ([^ rnn from jnst below Mf to Mew and the 


of the axial cnrrent 
operators in Eq. 


51 


rnn from jnst below Mew to Mqcd- Mnch of the literatnre concerns spe- 
cihc BSM models, where particnlar operators are singled ont. A model-independent analysis 
to one loop, focnsing on flavor-conserving interactions of the hrst generation, was carried ont 
in Ref. 


37 , and the more general sitnation is nnder stndy 52 . 
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For Mf ~ a few TeV, the dominant effects in the evolution of coefficients are due to the 
strong interaction, and can be written in terms of the numbers of colors and flavors rif, 
the Casimir Cp = — l)/2iVc, and the beta function for gg = -y/dvi^, — 2n/)/3. 

The EW running for some of the operators has also been considered in Ref. , where earlier 
references can be found. The collection C = {Ci,C 2 ,C 3 )'^ of d = 6 operators obeys the RG 
equation with a matrix 7 of anomalous dimensions, 


dCjfx) 

dln/i 




(y.o 




7dip 7mix 713 

7ps 0 
0 733 


( 11 ) 


Here Ci = (do, ds, do, d^, d^)'^ and C 2 = (ImSi, ImSg)^. The renormalization and mixing of 
qEDMs, qCEDMs and gCEDM is described by 


7dip = 


/ 8 Cp 
0 
0 
0 

V 0 


0 

SCp 

0 

0 

0 


—8(7 p 
0 

16Cp - 4iVe 
0 
0 


0 

—8 Cp 
0 

16Cp - 4iVc 
0 


while 


7ps = 2 


/3o - 2(3 + A/N,)Cp -2(1 + l/N,f(N, 
4(l + 2/iV,) (]^ + 2{Cp-l- 


0 \ 

0 

2 N, 

- 2 eNc 

/^o + + 2nf j 


-2)Cp 

2/Nl) 


( 12 ) 


7n 


-5 + 3£ 3-5£ -4 -4e 

(-5 + 3e)Cp (3-5e)Cp 2(N^-2Cp) 2e(N,-2Cp) 


(13) 

(14) 


contain, respectively, the anomalous dimensions of singlet and octet PS4QOs [^|^[^, and 
the mixings of dipoles with PS4QOs 37,55 . The hrst four columns of 7 dip are known to two 


loops 56 


Above Mew; G 3 = (g^ImEi, 0',ImY''^,IinY'‘^)^. The dipoles mix only with the gluon- 


Higgs operator 0', whose diagonal entry vanishes [^,57 , the right-handed current operator 
ImSi does not require renormalizatio n p^ , and the quark-Higgs couplings renormalize 

multiplicatively as the quark masses |37], so that 


7i3 — 


27r^v^ 


0 (0 0 1 e 0)'^ 0 0 


733 — —6Cp 


0 0 0 0 
0 0 10 
0 0 0 1 


(15) 


In these expressions we assumed that the qCEDM is proportional to the quark mass, and 
neglected terms with higher powers of quark masses. 

Below Mew, the quark-Higgs and gluon-Higgs operators disappear from the operator basis, 
and two new operators, the LR4QOs, appear: now C 3 = (ImSi, ImSg)^. These operators do 
not mix with the remaining operators and we have [^|54l[55] 


7i3 = 0 0 


733 — 2 


/3o 


/3o 


-3Cp/N, 
3N,(1 - 2/Nl) 


(16) 
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Besides the tree-level matching to coefficients of the Lagrangian ([^, other tree-level contri¬ 
butions exist. For example, the quark-Higgs operator contributes to chiral-symmetry-breaking 
four-quark operators. However, these contributions are suppressed by additional powers of the 
EW scale, and are effectively of higher dimension —we refer to Ref. 37 for a discussion. Most 
loop corrections involve powers of Oem, which for simplicity we ignore. The gCEDM receives 
threshold corrections from the top CEDM, and, at lower energy, from the bottom and charm 
CEDMs |53 


to a shift 


At one loop, these contributions are hnite, and, with our conventions, amount 

oisirng) 


6dw{mQ) = 


Stt 


-dgirng) {Q = t,c,b). 


(17) 


In the absence of the PQ mechanism, heavy quark CEDMs are also constrained by the large 

When the Higgs is integrated out, 9' generates an 


58 


radiative corrections they induce in 9 
0{as) threshold correction to the qCEDM, whose effect is however smaller than the qCEDM 
induced by running. Other important threshold corrections involve the top Yukawa coupling, 
and arise at two loops, through Barr-Zee type diagrams 


59 


The overall effect of the RG is to modify the coefficients by factors which are typically of 
0(1) (except for ImEi which can get enhanced by almost an order of magnitude [^). In Fig. 
we illustrate effects of RG evolution. We consider cases where at the scale Mf only the gCEDM 
(continuous lines) or the top GEDM (dashed lines) exists. For the sake of illustration, we took 
dt{Mf) = —lQQdw{,Mf). In the hrst case, we see that RG evolution reduces dw to about 20% 
of its original value, while generating a qGEDM and, to a lesser extent, a qEDM. While we 
plotted only the isoscalar components, also the isovector qGEDM and qEDM are generated, 
in the proportion d^/do = d^/do = e ~ 1/3. The tGEDM contributes to the gGEDM at the 


top threshold, Eq. (17), and generates light-quark qEDM and qGEDM through RG evolution. 


Although the induced gGEDM is a factor of 1000 smaller than the constraint from the 
neutron EDM is still about two orders of magnitude stronger than the direct bound from ti 
production at the LHG [^ . 

The main outcome of this analysis is that at /r ~ Mqcb the best organizational principle 
is given not by canonical dimension, but by the effective dimension inherited from the SM 
Lagrangian, Eq. ([^. Thus all the coefficients in Eq. ([^, including the q(G)EDMs, have sizes 
consistent with effective dimension six 


35 



, I w 





(18) 


where do, 3 , do, 3 , w, and ^ are eight real parameters encoding details of the BSM physics. 
Naive dimensional analysis (NDA) 1^,62 suggests that these parameters are 0(1). In con¬ 


trast, six other four-quark operators which are invariant under SU{3) x U{l)em 55 


originate from either higher-dimensional operators at Mew or higher-order electroweak effects, 
and are suppressed with respect to the PS4QOs and LR4QOs kept in Eq. 


35 37 
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Figure 1: Examples of the RG evolution of d = 6 operators: qCEDM do, qEDM do, and 
gCEDM dw (in units of as functions of the scale jj, (in GeV). The only nonvanishing 

operator at Mf ~ 1 TeV is taken to be the gGEDM (continuous lines) or the top GEDM 
(dashed lines). For visibility, the qGEDM, qEDM and gGEDM originating in the tGEDM are 
multiplied by 100. 

3 ^ at the Hadronic Level 

The^ interactions among quarks and gluons discussed in Section [^translate, at low energies, 
into ^ interactions involving the lightest mesons and baryons, which we sketch in this section. 


3.1 Chiral Symmetry and Low-Energy Interactions 


At a momentum comparable to the pion mass, Q ~ tiIt, Mqcd, the implications of the La- 
grangian in Eqs. (1^ and (1^ for the interactions among pions and nucleons are described by 


Ghiral Perturbation Theory (xPT) [^ 65 and its extension to arbitrary number of nucleons, 
Ghiral EFT |^-68 . The special role of the pion is a consequence of the invariance of Eq. 
Q (for e = 0) under the chiral symmetry SU{2)l x 517(2)^ ~ S'0(4), and its spontaneous 
breaking to the isospin subgroup SU{2)v ~ 50(3). Pions emerge as Goldstone bosons whose 
interactions are proportional to their momenta, which guarantees that low-momentum ob¬ 
servables can be computed in a perturbative expansion in powers of Q/Mqcd- Quark masses 
and other interactions explicitly break chiral symmetry, give pions masses, and induce non¬ 
derivative pion couplings, but the breaking is small and can be incorporated in the expansion. 
The chiral Lagrangian contains an infinite number of operators which can be grouped using 
a “chiral index” A = d + f /2 — 2 > 0 that increases with the number d of derivatives and 
chiral-symmetry breaking parameters, and the number / of fermions. Assuming NDA, the 
prediction of any observable at a given accuracy in Q/Mqcd requires the consideration of only 
a hnite number of operators, up to a certain A. 

Ghiral symmetry is realized nonlinearly in the chiral Lagrangian 69 71 , whose construc¬ 
tion via chiral-covariant objects is well known 72 . The choice of fields is arbitrary, and 
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for definiteness we employ a stereographic parametrization of the isospin-triplet tt, and an 
isospinor N = (pn)'^ for the nncleon. In this case the chiral-covariant derivatives take the 
formD^TTi = D~^{df,Sij + eAf,e3ij)7rj andV^N = [df,+iT-7TxD^TT/F^+ieAi,{l + T 3 )/ 2 ]N, with 
Fjr ~ 186 MeV the pion decay const ant0 and D = 1 + Isospin-invariant objects are 

antomatically chiral invariant; chiral-variant interactions are bnilt with the S'0(4) transforma¬ 
tion properties of the corresponding quark-glnon interactions. For example, electromagnetic 
interactions in Eq. Q break chiral symmetry as an antisymmetric S'(9(4) tensor [^. Be¬ 
cause mAT ~ Mqcd, Chiral EFT is well defined only for non-relativistic nucleons. The Q/rriN 
expansion can be made consistent with the Q/Mqqd expansion with heavy nucleon helds [^ , 
when the nucleon mass m^v is removed from propagators and the Dirac structure simplihes 
to the nucleon velocity and spin S^. Resummations of the Q/rriN expansion are currently 
popular but they do not decrease the overall theoretical error. Chiral EFT can be extended 


to the Delta-isobar region using explicit helds for the Delta 75,76 and the Roper 77 , but 


most ^ applications have so far been limited to nucleons. Chiral EFT reproduces the (B)SM 
S matrix because it includes all operators consistent with RG invariance and the symmetries 
of quark/gluon interactions at /r ~ Mqcd- 

Equation (|^ is represented at lowest chiral index as 


p(o) _ 

JL^rp - 


■D>^7z + N{zvV 
FCsNNNN + CvNtN ■ NtN, 


‘^S^T ■D^tz\N 


(19) 


where qa = 0{1) ~ 1.27 is the pion-nucleon axial coupling and Csy are LECs related to the 
two nucleon-nucleon scattering lengths. Higher-index interactions are constructed with further 
derivatives and nucleon helds. We assign a chiral index 3 to aem/47r, powers of which purely 
hadronic operators generated by the integration of hard photons — e.g., the electromagnetic 
pion mass splitting = (!7(Q!emd7QCD/47r) ~ 1260 MeV^ — are proportional to. 

The properties under chiral symmetry of the operators in Eqs. (§ and (|^ dictate how 
to incorporate them in the chiral Lagrangian, and determine the relative importance of ^ 
couplings 33 35 . The average quark-mass term in Eq. 0 breaks chiral symmetry as the 


fourth component of an vector, and induces interactions proportional to powers of 

fhr{6). The quark-mass splitting and the 6 term transform as diherent components of a single 
other ^(^(d) vector, which implies that their hadronic matrix elements are directly related 
They generate interactions proportional to powers of rher~^{d) and rhr~^{6){l — 


e^) sin 0/2, respectively. For momenta Q ~ m^, and taking er~‘^{6) = 0{1), quark-mass terms 
are paired with chiral-symmetric operators by making d count powers of the pion mass as well. 
The lowest terms stemming from Eq. 0 are 




277^ 

FfD 


NN+^^N 


Ts — ^ ^ 4- 4—^ sin 9 ) N. 

F^D \F^ 2e J 

( 20 ) 

The pion mass is = 0{MQCj:,mr{9)), the nucleon sigma term is Am^f = 0{m^/M qqj)), the 
neutron-proton mass splitting from the quark-mass difference is = 0{er~‘^{6)m‘^/MQCT)), 


7(3 1 _ 

F^ 2e 




sin0 


^ Note that the pion decay constant is frequently defined in the literature as F = Fr/2 or F = 
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and pion-nucleon interactions —PT charge-symmetry breaking (CSB) for an even nnmber 
of pions, for an odd number— are determined by the matrix element that enters 6m n. 
Similar relations exist at higher orders —for example, involving the hadronic contribution 

to the pion mass splitting 6ml = 0{6m %)— but the link between PT CSB and 9 operators 
quickly becomes superfluous. 

The d = 6 operators in Eq. ([^ have different transformation properties under chiral 
symmetry, and the connections with PT operators are more tenuous. The qCEDM and 
qEDM break chiral symmetry as 50(4) vectors and induce both isospin-conserving and isospin¬ 
breaking interactions, which are in general of the same size because after aligning the 9 
term there is no longer freedom to eliminate the qCEDM and qEDM isovector components. 
The important difference between qCEDM and qEDM is the suppression for the latter of purely 
hadronic operators by powers of aem/47r. The LR4QOs also break chiral symmetry (and isospin 
in particular), but as 50(4) tensors, and the relative importance of their interactions is similar, 
but not identical, to the isovector qCEDM. The gCEDM and PS4QOs do not break chiral 
symmetry and cannot be distinguished purely on the basis of their symmetry properties, more 
information about their matrix elements being required. For chiral-invariant operators, chiral- 
breaking interactions, like non-derivative pion-nucleon couplings, are suppressed by factors of 
the quark masses. 

The PT pion, pion-nucleon and multinucleon chiral Lagrangians are known up to A = 

The Lagrangian from 9 and d = 6 operators was built in great 
For nucleon and light-nuclear EDMs at LO it is sufficient to 
consider a subset of the interactions discussed in Refs. ISSHSSl: 


4, 3, 3, respectively 

78 

00 

detail in Refs. 33 

36 

82 

00 


Cf = -^N {goT-TT + gi7r3)N-2N (do + diTs) S^^Nv'^F^ 

■^71 

A 




-^TTgTT^ + CiNNd^ (NS^^N) + C2NTN ■ (NSf^rN) . 


( 21 ) 


The operators that couple a neutral pion to the vacuum (pion tadpoles) can be eliminated 


order by order in favor of the interactions remaining in Eq. (21). Each of the operators in 


this equation has chiral partners, which we do not display explicitly. The interactions in the 
second line are only needed at LO for LR4QOs, PS4QOs, and gCEDM, while for qEDM only 
EDM-type operators in the first line are important. The coupling constants ^o,i) do^i, Ci^ 2 , 
and A are discussed in the next section. 


3.2 Sizes of T Couplings 


The LECs in the ^ hadronic Lagrangian, Eq. (21), are (approximately) linear functions of 
9 and of the coefficients of d = 6 operators in Eq. ([^. We impose the constraints of chiral 
symmetry and use NDA for an estimate of the scaling of the LECs. In order to connect high- 
energy observables to EDMs, it is crucial that the coefficients of proportionality be known 
accurately. Going beyond NDA requires nonperturbative techniques, and we review here 
information that can be extracted using symmetry, LQCD simulations, and QCD sum rules. 
(A more comprehensive review of calculations of these matrix elements can be found in, for 
example. Ref. (4^.) 
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The LECs Qq and gi are isoscalar and isovector ^ pion-nucleon conplings, which induce 
important long-range contributions to the nucleon EDM (Section and to the nucleon- 
nucleon potential (Section [6^. One expects 33,35 


% = o 


m; 


^QCD 


e, mldo, emlds, mldw, Mqcd , ( 22 ) 


4 4 

eml - em: 


^ ^ \ I 1 Mqcd 1 • (23) 

QCD QCD 

Because ^o,i Eire LECs of chiral-breaking interactions, for most sources the quark masses 
appear. The contribution of isoscalar operators to gi is suppressed by factors of rri^/ 
but as we discuss below gi can still be important in systems, like the deuteron, for which the 
contribution of g^ vanishes. Isovector operators generate both go and gi at the same order, 
even though is affected by the quark mass difference. For all the operators in Eq. (|^, a 
third non-derivative pion-nucleon coupling, (Jj 2 /FT,)NTi^T‘iN, is suppressed, and contributes at 
the same level as derivatives^ couplings which we neglect. 

Important consequences of chiral symmetry and vacuum alignment are the survival of the 
three-pion interaction A and modihcations of So,i proportional to it. Both So,i receive tree-level 
contributions. Moreover, at one loop A gives a signihcant contribution to the pion-nucleon 


FF. Although formally an NLO effect, the loop is enhanced over NDA by a factor of dvr 35 


It endows the FF with a certain momentum dependence, whose effects on nucleon EDMs have 
not been studied. For light nuclei, they were found to be small [^. We can capture the 
momentum-independent effects by redehning Si- Thus, 


9o — 9o + 


6m 


N 


mt 


A + 


S'! — S'! + 2 


/ ArriAf _ 15 ^^ 771 ^ 
\ m‘1 


A + 


(24) 


where are the couplings before alignment. The main remaining, explicit A contribution is 
a tree-level three-nucleon potential (Section |6.l[ ). Only for LR4QOs is A an LO effect, while 
gQ is higher order. In this case, we can eliminate [35], 


A = 


m: 


6m 


9o — ^ "^qcd) 

N 


(25) 


The pion-nucleon couplings have been best studied when links to PT quantities through 
chiral symmetry are useful: 


For 9, a comparison between Eqs. (20) and (21), and similarly for higher-order terms, 
yields ^ 33 , 


Ah 


F^ 

El 

F^ 


6mN 1 — 


Fn 2e 
Atutv lbg\m. 


sin 6 = (15 ± 2) ■ 10 ^ sin S, 
6mi 1 — 


(26) 


mt 


IGttFJ 


sinS =-(4±3) ■ 10"^sinS, (27) 
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using the values from LQCD, e = 0.37 ± 0.03 
and AmAT = —63 ± 9 MeV 
MeV^ 


92 


5mN = 2.39 ± 0.21 MeV 89,90 


91 , and from xPT fitted to meson data, 5m^ = 87 ± 55 


This value for qq is somewhat smaller than the NDA estimate. The relation 
between and 5m]s[ is violated by terms of 0{m ^/an effect of the same size as 


the uncertainty in Eq. (26). gi is related to PT operators that contribute not to baryon 
masses, but to pion-nucleon scattering observables —however, only at an order beyond 
the current most precise analysis 


93 


For qCEDM, analogously 35,94 


go = SruN 


druN do Ami do 
6mN Co ml cq 


gi = 2 Am 


Am 


N 


N 


Am 


N 



where co ,3 are the coefficients of chromomagnetic operators analogous to the chromo¬ 
electric operators in Eq. ([^, Ami, Am^, and 6m n are corrections to the pion mass, 
nucleon sigma term, and nucleon mass splitting due to Cq, Cq, and 63 , respectively. The 
evaluation of the corresponding matrix elements is currently being pursued by lattice 
collaborations 


95 


rules. It is found 


. At the moment, the best estimates on go^i come from QCD sum 
11 that gi is larger than go by a factor of about 5, and 


= -{20tTo) ■ 10 -=" {27iF^fdo. 

r-jT 


(29) 


The large error is due to cancellations in Eq. (28), which complicate estimates of these 


couplings. In terms of dimensionless quantities, the numerical factors in Eqs. (26) and 


(29) are not very different, as one would expect from Eqs. (22) and (23). 


The parameters hop represent short-range contributions to the nucleon EDM (Section 
and have expected sizes 33,35 


do,i = O 


m: 


^QCD 


6, m^di, m^di, Mq^-^dw, Afg^^ImSa, Mq^plm,; 


M< 


QCD 


(30) 


Chiral breaking intrinsic to the electromagnetic interaction ensures that no extra factors of 
quark masses are needed beyond those appearing explicitly in Eqs. (|^ and ([^. Because some 
of this breaking involves isospin, all sources induce isoscalar and isovector components of the 
same size, do.i are not hxed by symmetry but can be extracted from LQCD calculations of 
the nucleon EDMs, as we discuss in Sect. 


5.2 


The remaining couplings in Eq. (21), Ci^ 2 : are less well known, and have not been studied 
extensively. They represent short-distance contributions to the nucleon-nucleon potential and 
their main phenomenological impact is through light-nuclear EDMs (Section]^. They are 
chiral-invariant interactions that appear without suppression, and at LO, only for 


chiral-invariant operators 35 


Cl ,2 — O I ((ii4/,ImSa) 


M, 


QCD 


F 2 


(31) 
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For this reason, contact terms with different isospin structures are not needed at LO. 


Equations (22 )-(31) show how the different properties of 6 and d = Q operators under chiral 


symmetry lead to different hierarchies among the couplings in Eq. (21), which has profound 


consequences for the EDM of the nucleon and light nuclei. Taking into account Eq. (25), these 
observables are expected to be described in LO in terms of the six independent LECs go,i, do,i, 
and Ci^ 2 - In Sectionj^we find that the nucleon EDM depends mostly on ^o,i and do,i) whereas 
in Section we show that the six LECs enter nuclear EDMs. By contrast, the traditional 
description of nuclear ^ is largely based on the three non-derivative pion couplings 
and g 2 , which is sub leading for all d = 6 sources. 


96 


9o,i 


4 Peccei-Quinn Mechanism 

So far, we have assumed that there is no particular reason behind the small value of 6. A very 
elegant way to obtain a small 9 dynamically is the PQ mechanism (^, which is discussed 
in EFT in Ref. 97 . The approximate, spontaneously broken U{1)pq symmetry is realized 


nonlinearly, with SM fields invariant and the pseudo-Goldstone boson, the axion a 18,19 


changing by an additive constant, a —?• a + c. The symmetry is explicit ly b roken by the 


anomalous coupling to GG, which can be eliminated, as we did in Section 2 T, in favor of a 
complex quark mass term —except that now the axial rotation depends on a. 

After vacuum alignment, which in this context is equivalent to the diagonalization of the 
pion-axion mass term, the axion Lagrangian reads 


jC-ax = ^df,ad>^a + [Go + {Gi + e r ‘^{9 +a/fa )) Ts] 757 ^^ d>^a 


(32) 


-mr{9 + a/fa) qq + emr ^{9 + a/fa) q 


(1 — — 

t -3 + sin(6' a/fa) ^75 


q + ... 


where fa is the axion decay constant and the two couplings Coq are model dependent. When 
chiral symmetry is broken, the term proportional to qq generates the axion potential V{9 + 
a/ fa), whose minimum (curvature) determines the axion vev (a) (mass). If 9 is the only source 
of the potential is minimized hj 9 + { 0 )/fa = 0. The presence (in the “... ”) of higher¬ 
dimensional QP operators that break chiral symmetry affects this potential; for example for 
the qCEDM, 


V{x) = 


m: 


rfx) 


1 - 


Ami 

2coml 


do -|- eds ) r ^{x) sin a; 


Minimization results in an induced angle 


0ind = 0 + ^ = - 
Ja ^ 


2 Ami 


£2 coml 


do £ d^ 


(33) 


(34) 


of 0{M^q^/ similar to other d = 6 operators. The consequence to low-energy dynamics 
is that the coupling go induced by the qCEDM receives another correction. 


(5'o)pq — 90 + SmN - 


2e 


9\nA — 


6m 


N 


ec3 


S + 


C 3 


Ami 


Co 


do 


(35) 
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which cancels the contribution from d^. Similar relations can be worked out for LR4QOs, 
resulting in a vanishing Qq. 


5 Nucleon Electric Dipole Moment 


The study of the nucleon EDM in yPT has a long history starting with Ref. [^, where 
the leading pion-loop contribution to the neutron EDM induced by 6 was computed. The 
calculation was later extended to the radius 98 and then full 99 electric dipole form factor 


(EDFF), and to NLO 100-|102 


computed to NLO in Refs. 35 103 


the EDFF is interesting. Due to Schiff’s theorem 104 


The nucleon EDFF generated by the d = 6 operators was 
In addition to the EDM, the momentum dependence of 
the nucleon EDM does not contribute 


to the EDM of atoms in the nonrelativistic limit, and the first non-vanishing contribution is 
induced by the EDFF radius 

to guide the extrapolation of LQCD results in both pion mass 


Furthermore, the EDFF presented in Section |5.1| can be used 

107| and momentum. 


105 


as 


described in Section 5^ The role of strangeness is examined in Section |5.3 


5.1 Chiral EFT 

yPT allows for the calculation of low-energy observables in a controlled perturbative expansion 

A review, including the 
The toroidal 


63 


in powers of Q/Mqcb, where each loop contributes 
PT electric charge and magnetic dipole FFs, can be found in Ref. 108 


can be written as 99,103 


dipole FF can be found in Refs. 109, 110| . The component of the electromagnetic current 


J^{q,K) = 2{Fo{Q^) + Fi{Q^)Ts) 


S>^vq-S-qvf^ + — (S^q -K-S- qK^) + ... 
niN 


, (36) 


where q = p — p' and K = (p -|- p')/2 in terms of the nucleon momentum in the initial (final) 
state p {p'), and > 0. Here Fq{Q^) {Fi{Q‘^)) is the isoscalar (isovector) EDFF of 

the nucleon, 

F,{Q^) = d,-S'Q^ + H,{Q^), (37) 


where di is the EDM, S'' the Schiff moment, and Hi{Q‘^) accounts for the remaining de¬ 
pendence. 


dn T dj, 


— dn -|- 


32 

35 

100 

103 

e9A9o 

- -— TT 

3m^ 


4m AT 


1 + #^ 1 - 


6m 


N 


dp d^ ^9 a9o 

— = di + 


{27rF^ 


Spo 

^ dvr m 


m^ 


L-ln^ + 

/i^ 4 mjsf 


1 + 1^ 


6ml 


mt 


(38) 

(39) 


where L = 2/(4 — d) — 7 ^ -|- log 47 r, with d the spacetime dimension and 7 ^ = 0.557... the 
Euler constant. For illustration, the leading diagrams are shown in Fig. The nucleon 
magnetic dipole moment (MDM) couples to the electric field in the A = 2 Lagrangian, an 
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Figure 2: Leading diagrams for the nucleon EDM. Solid, dashed and wavy lines represent prop¬ 
agation of nucleons, pions and photons, respectively. Dots (squares) denote T interactions. 
For simplicity, only one possible ordering is shown. 


effect oc rrij/ 


It does not contribute to the EDM before next-to-next-to-leading order (N^LO), 
contrary to uncontrolled calculations based on a “relativistic” chiral Lagrangian [m], The 


subset of 0{l/mj^) corrections that are proportional to the MDM is small 112 


From the NDA estimates of Section 3.2[ we see that the relative importance of the terms 
in Eqs. (38) and (39) depends on the source: 


Fig. I^a), dominates, while loop contributions. Fig. [^b-d), appear at least at N^LO 
for qEDM, gCEDM, and PS4QOs. The full set of N^LO corrections was considered in 
Ref. 


103 . 


go appears at LO for sources that break chiral symmetry and do not contain photon 
helds: 6, qCEDM and LR4QOs. Its one-loop contribution is purely isovector, and the 
dependence on L and fi is eliminated solely by di(/i). Non-analytic corrections to the 
isoscalar EDM are finite and suppressed by Despite the extra factor of vr with 

respect to NDA, they are about 10% of the leading loop. However, one cannot exclude. 


on the basis of Eq. (30), that do is present at the same order. Without further dynamical 


information, one cannot state, as sometimes in the literature, that the nucleon EDM is 
isovector at this order. 


• Qi contributes through recoil corrections in the nucleon propagator and pion-nucleon 
axial couplings, and only affects the proton EDM. In the case of isoscalar sources, gi is 
formally suppressed by Even for 6, where go is a factor of ten smaller than 

expected by power counting, gi gives a correction to the proton EDM which is only a few 
percent of the leading loop. In the case of isospin-breaking sources, ^o,i appear at LO, 
but ^0 arises from pion tadpoles and is proportional to Sruisr, which makes it numerically 
smaller than gi [^. For the isovector qCEDM and LR4QO, then, the nucleon EDM is 
likely to receive its largest numerical contribution from the LECs do,i, with non-analytic 
terms entering at the 30% level. 


For concreteness, we specify Eqs. (38) and (39) for 9, where go is well determined by 
Eq. (26) and gi gives only an N^LO effect. Setting the renormalization scale fi = rriN and 
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neglecting the nnmerically small isospin-breaking contribntions, 


dn - {do - di) = 

dp ~ {do + = 


egAdo 

“(27rF^)2 

egAh 


In 


m 


N 


mi 


71 m-n 

+ 


2 rriM 


{27rF^ 


In 


m 


N 


mt 


+ 277- 




rriN 


9 = -{1.99 + 0.12) ■ 10-3 sin 9 e fm, (40) 
0 = (1.99 + 0.46) ■ 10-3 sin 0 efm, (41) 


showing that, especially for the nentron, convergence of the SU{2) chiral expansion is good. 
Assnming that there are no hne-tnned cancellations 32 between doA, which are analytic in 
ml, and non-analytic contribntions, the current bound on dn allows to p ut a bound on 9. 
Similarly, using the NDA expressions for ^o,i and do.i given in Section 3.2 and assuming no 


cancellations, one can bound the coefficients of d = 6 operators 103 


dj MqQYydi, MQQ-^di < 10 


-10 


^QCDdw, 




-12 


QCD-* 


(42) 


The weaker bound on 9, qCEDM and qEDM reflects the proportionality to light-quark masses. 


If the dimensionless^ parameters in Eq. (18) are 0{1), the bounds in Eq. (42) suggest new 
physics at a scale Mr > 100 TeV. Once NDA can be replaced by more accurate determinations 
of matrix elements in LQCD, these bounds can be reliably translated into bounds on the 
dimensionless parameters at Mf using the RG of Section 


2.2 


as sketched in Ref. 37 


While certainly an exciting evidence of new physics, the observation of neutron and proton 
EDMs would not, with the current theoretical status, be sufficient to clearly identify the source 
of new physics, and in particular to disentangle the effect of 9 from higher-dimension operators. 


Without further input from nonperturbative techniques, Eqs. (40) and (41) show that, for all 
operators in Eq. (1^, dp^n contain at least two unknown LECs. An interesting independent 


observable is the momentum dependence of the EDFF. To NLO 35,98,99,102,103 , the square 
radii are 


, ^ TTegAgpSmN ^ egAgo 

° 12(27rF^)2m3’ ^ 6(27rF^ 


|2^2 


1 - 


Svr 

4 mN 


Sml 


mt 


(43) 


For 9, qCEDM and LR4QO, radii arise at the same order as the EDM, are finite in the chiral 
limit and approximately isovector. They are dominated by contributions at the scale m^^, 


which are known except for the ^ parameters. For example, for 9 102 

-5.0 ■ lO-®sin0 efm3, 




S[ = 6.8 ■ 10 3 sin 6* efm' 


(44) 


In contrast, for qEDM, gCEDM and PS4QOs, radii arise at N^LO and scale as with 


respect to the EDM 103 . The functions Hi{Q^) from Eq. (37) can be found, to NLO, in 
Refs. 35,99,102,103 . Although an experimental measurement of the momentum dependence 
of the EDFF is not going to happen any time soon, the extrapolation —)■ 0 of LQCD EDFFs 
for 9 and qCEDM would allow to extract go and do.i at the same time. 


5.2 Interplay with Lattice QCD 

The best tool to determine the LECs do,i that contribute to the nucleon EDM is LQCD. 
Unfortunately there are virtually no LQCD results for the d = 6 sources, most work having 
focused on 9. A collection of early results can be found in Ref. 
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A recent LQCD evaluation of dpn from 6 can be found in Ref. 


114 


Simulations were 

performed at two values of the pion mass, = 330,420 MeV, with domain-wall fermions to 
suppress the lattice artifact of chiral symmetry violation. Results, which still have a significant 
statistical error, are extrapolated linearly to Q = 0, where they give EDMs compatible with 
zero. Ref. 107 improved the finite-volume corrections of Refs. 105,106 and used these LQCD 


results to extract dop with SU{3) yPT at NLO. We have fitted the same EDM points but 
with the SU{2) formulas, Eqs. (38) and (39), and neglecting finite-volume corrections. We 
hnd good hts, as exemplihed for the isovector component in Fig. where the best £t results 
and Icr uncertainty obtained from Eq. (39), with go hxed to the value ( |26| ), are compared with 
results of fitting di only and neglecting the chiral loop, that is, setting go to zero. Clearly, 
lattice EDM data show an essentially linear dependence on with no clear sign of the chiral 
log. For the short-range LECs at the scale g, = 939 MeV, we hnd, factoring out the pion-mass 
dependence. 


do{g) = (-0.04 ±0.45) 


ml 


{27rF^ 


esind, di(/i) = (0.05 ± 0.45) 


m: 


{27rF^ 


esind. 


(45) 


which is not far from the NDA estimate of Section 3±, once one considers the large uncertain¬ 
ties. Extrapolating to the physical pion mass, again using Eq. (26), 

d„ = —(2.4 ± 1.5) ■ 10“^sindefm, dp = (2.5 ± 1.5) ■ 10“^sindefm, (46) 

which are consistent with the more sophisticated analysis of Ref. |107| . After this review was 
completed, an LQCD calculation with clover fermions and imaginary 6 has appeared 115 , 
which gives a similar value for dn but with smaller error bars . 

Eventually, the extrapolation in could also be made with the full yPT EDFF, Eq. 



The uncertainty on the LQCD evaluation of doq is quite large and compatible with zero. 
Nonetheless, for 9 LQCD is now competitive with other nonperturbative techniques. For 
comparison, from QCD sum rules [117 - 119] , 

|d„| = (2.5± 1.3) ■ 10-¥efm. (47) 

QCD sum rules also give the best available estimate for the qCEDM [120] , 

\dn\ = (1± 0.5) (27rF^)2|0.9do-0.3 41 -lO-^efm. (48) 

Among d = 6 sources, the qEDM is particularly simple, since to a good approximation 
dp^n arise solely from do,i with small isospin-breaking corrections. From Eqs. (|^ and ( [2l| ), the 
photon couples in LO to the isoscalar and isovector tensor charges dehned by {N\qap_i,q\N) = 
‘ig'^e^’^^'^VpNSa-N, {N\qap^Toq\N) = Agfe^’^^'^VpNSfjToN, resulting in 


dji ± d^ 


dp dji 


dome 


= gjme 


do , 

- dri 

3 

do± y ) = (0.51 ±0.04) me 


= (0.29 ±0.03) me ( y + ^^3 


j da 

do + J 


(49) 


(50) 
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Figure 3: Isovector component di of the nucleon EDM (in units of sin^e fm) as function of 
the pion mass squared (in units of GeV^). The black lines and gray shaded area show the 


best £t and la uncertainty obtained by htting Eq. (39) to the LQCD results of Ref. 114 
represented by the red points with error bars. The blue lines and blue shaded area are the 
corresponding £t with qq = 0. The dashed vertical line indicates the physical pion mass. 


where we used the recent LQCD results 121-123 for the proton tensor charges at p = 2 GeV, 


and quark mass and qEDM should also be evaluated at this scale. These numbers are in good 
agreement with the NDA expectation. Existing extractions of the proton tensor charges from 
experiment 124,125 are at too high a for yPT. Some of the issues affecting the lattice 


implementation of higher-dimensional operators are discussed in Ref. 126 


Our simplified analysis of LQCD data, which omitted correlations, finite-volume effects 
and systematic errors, suggests that there can be a rich interplay between LQCD and yEFT. 
It will be interesting to see whether the substantial reduction of the LQCD uncertainty to the 


10% level expected in the near future for 9 114 will show a clear signal of an EDM and the 


presence of a chiral logarithm, allowing an extraction of ^o- More generally, LQCD calculations 
of matrix elements for other sources would £11 a gaping hole in the road from BSM physics to 
EDM experiments. 


5.3 Role of Strangeness 

Much of the low-energy ^ literature includes the strange quark s explicitly. Yet, the inter¬ 
mediate value of its mass poses a problem: it prevents integrating s out at a perturbative 
scale, as it is done for the t, b and, to a lesser extent, c quarks, but leads to poor convergence of 
the yPT expansion via a relatively large kaon mass mx 127 . Although the best way to deal 
with s (and maybe also c) is through nonperturbative techniques, symmetry considerations 
also give some insight. 

9 e£ects in SU{3)l x SU{3)r {U{3)l x U{3)r if one wants an explicit connection to the 
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U{1)a anomaly) xPT 82 86 are proportional to 


_ mdruums __ 

ms{md + rriy)+murnd 2 ^ 


m 


1 + —(l-£^) + 

m.o 


( 51 ) 


Althongli go and Smisi receive, starting at order O^mK/Mqci)), large corrections from kaon 
and eta meson loops (which show little sign of convergence), these corrections affect go and 
6mN in exactly the same way np to N^LO and are already acconnted for throngh Eq. (26) 


when extracting 6m]sr from LQCD simnlations with dynamical strange qnarks. At N^LO, in 
both SU{2) and SU{3) xPT there are short-distance contribntions to go that do not affect 
6mN, bnt, since they are not proportional to m^, their nnmerical effect shonld not be larger 
than the nncertainty qnoted in Eq. (26). The isospin-violating conpling gi receives a tree- 


level contribntion proportional to the n-g mixing angle, which, despite being formally LO, is 
snppressed with respect to go by fh/mg and is nnmerically of the same size as the estimate of 


Eq. (27). The hrst contribntion not snppressed by powers of m/mg appears at N^LO, as in 


SU{2) yPT. Thns, strangeness shonld not signihcantly affect the valnes of the pion-nncleon 
conplings from 6. 

The nncleon EDM has been compnted in this framework to varions degrees of sophistica- 
At one loop, there are additional nncleon-kaon and nncleon-eta conplings. 


32,101,106 


tion 

which at LO are hxed by combinations of baryon masses, for example gN'ZK is related to the 
mass difference between the nncleon and the S baryon, mjv —ms. As for go, these relations are 
violated at N^LO. However, now the corrections are proportional to mg, and are nnmerically 
more important 128 . For the nncleon, pion and kaon loops give contribntions of approxi¬ 
mately the same size. However, NLO corrections to the kaon contribntions are as large as LO, 
casting some donbts on the convergence of the 517(3) expansion. On the contrary, the 517(2) 
expansion of the nncleon EDM converges fairly well. 

In the 7 = 6 sector, the basis in Eq. ([^ needs to be enlarged to inclnde s. Fonr additional 
qnark bilinears can be constrncted, the strange EDM (sEDM) and CEDM (sCEDM), and 
two A5 = 1 flavor-changing nentral cnrrents of the same form as the qEDM and qCEDM 
bnt containing a 7 and an s. There are also several more fonr-qnark operators, listed in 
Ref. 54 , which, however, does not fnlly acconnt for the constraint of SM gange symmetry. 


If one considers only AS = 0 operators, in addition to the fonr operators already dehned 
in Eq. (|^, eight fonr-qnark operators with two strange qnark helds receive non-vanishing 
matching coefficients at tree level, so that only half of the operators dehned Ref. 54 needs to 
be considered. The matrix elements of these operators were stndied (54) in the factorization 


approximation and in the qnark model, bnt more rigorons techniqnes are needed. 

The sEDM and sCEDM have received more attention. The effect of the sEDM dg on 
dp^n can be obtained from {N\sa^‘'s\N) = 2g'^e^^^^VpNSaN on the lattice: factoring ont the 
strange qnark mass and charge, as in Eq. dn + dp = —2gJmgds/3. gj involves sea qnarks, 
and is nnmerically more challenging to calcnlate on the lattice than g'^^. A recent, state- 


of-the-art LQCD calculation 121 Ends gj = 0.002 ± 0.011. With these large uncertainties, 
it is not possible to exclude that, due to the enhancement mg/md ~ 20, the sEDM gives a 
contribution to dn + dp of the same size as Eq. 


been addressed in 5f/(3) xPT and QCD sum rules In 


T 

re contribution of 

ffie sCEDM dg has 

111 

55 

1129). Refs. 

55 

129 

found that dg 


could give a large, possibly dominant, contribution to the nucleon EDM. In the presence of 
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PQ symmetry, dg does not contribute at LO to pion-nucleon ^ couplings, but affects g^EK 
and thus induces long-range contributions. The numerical value of the loop is small, but not 
so small as to compensate the enhancement due to the appearance of rUg in the coefficient of 
the sCEDM. It would be interesting to test the robustness of this prediction by going beyond 
LO. In general, the same issue as for 6 remains: the convergence of SU{3) yPT. 


6 T Moments of Light Nuclei 


The observation of neutron and proton EDMs will not provide enough information to clearly 
identify the leading^ source(s). ^ in light nuclei, being sensitive to different combinations 
of the couplings in Eq. ( [2l| ), is highly complementary, and it can be computed in the same 
theoretical framework as the nucleon EDFF with reliable accuracy. For these reasons, and with 
the additional motivation of the exciting experimental developments that might allow a direct 
measurement of EDMs of charged ions, the study of the EDMs of deuteron, helion and triton 
has received a lot of attention in the last few years, both from a phenomenological perspective 
and with Chiral EFT 


130-134 


87,135-137 . Very recently, the hrst model calculation of the 


®Li EDM has appeared 138 . For the deuteron, as seen in Tab. two other moments, MQM 


and TQM, are sensitive to^ and have been calculated for d <6^ sources in Refs. 135,139 
and |140 , respectively, but prospects for their detection are remote at best. ^ in neutron- 


proton [141 and neutron-deuteron 142 scattering have been revisited recently, but seem 


equally difficult to measure. 

is a small effect, and can be treated perturbatively on top of the strong nuclear interac¬ 
tions. The elements that enter EFT calculations in nuclei are discussed in Section 6.1[ while 
results for nucleon number A = 2,3 are reviewed in Section 6.2 The alpha particle {A = 4) 


has no^ moments, but T EEs of nuclei with a few more nucleons could be calculated within 
this approach. 


6.1 Nuclear Potential and Currents 

The T electromagnetic current can be written in the form of Eq. 
and in a straightforward generalization that includes MQM |l35 


for a spin-1/2 nucleus, 
and TQM 140| for spin-1. It 


receives several contributions, schematically illustrated in Fig. ^ The nuclear wavefunction, 
denoted by the shaded triangles, as well as the iteration of the PT potential Vp^, represented 
by the shaded blob, can be obtained by solving the Schrodinger equation. T', denoted by 
black squares, either affects the electromagnetic current (Fig. |^a)) or perturbs the nuclear 
wavefunction with an insertion of the^ potential (Fig. |^b)). 

The power counting of Chiral EFT systematically organizes the various contributions 


35 


136 


While in the one-nucleon sector a loop brings in a suppression of 


m 


problems with two or more nucleons a relative enhancement of ATmiM/Q compensates the 
loop suppression in diagrams that involve intermediate states with nonrelativistic nucleons 
only, where energies are 0{Q‘^/mt^) rather than 0{Q) as in yPT [^. The one-body contri¬ 
bution in Fig. I^a) represents the nuclear FF arising from the EDFF of the constituents, Eq. 
(37). The PT one-nucleon current contribution to Fig. |^b) contains an additional loop with 
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Figure 4: Contribution to nuclear EDMs. The large triangles denote the nuclear PT wavefunc- 
tion; the oval, iterations of the PT potential; the oval with an attached photon, a few-nucleon 
current; the black square, the^potential; and the large black square with an attached photon, 
a^ few-nucleon current. Other notation as in Fig. 


respect to the one-body contribution in Fig. |^a). These contributions scale as 

J. = 0(F„_,Q), A ^ O (^Q. eC^aF^Q, O . (52) 

\ TT TT -^QCD / 

Depending on the sizes of ^o,i compared to dn,p, insertions of two-body currents in Fig. [^b) 
can be as big as the one-body contribution to Fig. Ba) 

The various ingredients in Fig. |^can be obtained explicitly in terms of the ^ LECs. For 
EDMs, the most important ^ current is the sum of one-body contributions. 




dfi dp — 1 “ (^dp 


dn)r. 


(i) 


a 


(i) 


(53) 


where is the spin (isospin) of the nucleon i that interacts with the photon, and 

q is the (outgoing) photon momentum. Pion-exchange T' currents induced by go, which by 
power counting could be important for the deuteron EDM from 9, are given in Ref. 136 , but 
their contribution is found to be numerically small. 

The tree-level potential induced by the couplings in Eq. (21) has a two-body component 


34 , which is generated by one-pion exchange (OPE) with couplings go^i and by the short-range 


interactions Ci^ 2 , 


, 9Agi 


-FlC^ + 


“^gAgo 

kf + ml 


FtCo 1 


cr 


(b 


a 


(i)) 


,, , ,, 4’ + P) (<?® - a«) + ri‘' - r«> (J® + a®) 

and a three-body component [^, which is generated by the three-pion coupling A, 

k>j>i 

■ ki ■ kj ■ kk 
X 


{k^ + ml){k^ + ml){kl + ml) ’ 


(54) 


(55) 
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Table 2: Expected orders of magnitude for the neutron EDM (in units of e/Mgco), and for the 
EDM ratios proton to neutron, deuteron to neutron, helion to neutron, and triton to helion, 
for 9 and d = 6 sources. Q represents the low-energy scales m^, and with B the 

binding energy. (Adapted from Ref. 136|.) 


where ki is the momentum transferred from nucleon i 
couplings (7o,i is well known 


96 


The OPE potential generated by the 
It should be important for 9, qCEDM, gCEDM, PS4QOs, 
and LR4QOs, but for the 9 only qq appears at LO. The short-range terms, which can be 
thought as the long-wavelength effect of ^ heavy meson exchange, should be kept at this 
order for gCEDM and PS4QOs. The three-body force, which is new, appears at LO only 
for LR4QOs. The LO potential induced by the qEDM is given in Ref. 


34 , but it is of little 


phenomenological relevance. Chiral EFT provides a framework to go beyond LO by including 
two-pion exchange potentials, relativistic corrections, and OPE potentials from subleading 
couplings. For 9, N^LO corrections to Eq. (54) were discussed in Ref. 34 . The calculation can 


be immediately extended to the qCEDM, with the observation that TPE diagrams involving gi 
In the case of the LR4QO, the coupling A gives a relatively large NLO correction 

the largest part of which we absorbed in gi 


vanish 137 


to the isospin-breaking OPE potential 
(Eq. 


(24)). 


35,87 


6.2 Moments 


The estimate of Eq. (52), combined with the relative sizes of g^^i, Ci ,2 and A discussed 


in Section 3.2 lead to the expectation that for qCEDM and LR4QOs EDMs of light nuclei 


are dominated by the ^ OPE potential. For the gCEDM and PS4QOs, one-body, OPE and 
Ci ^2 contributions should be approximately equal, while if the qEDM is the only ^ source 
light-nuclear EDMs are well approximated by the EDMs of the constituents. For 9, the 
situation is more complicated. For most nuclei, OPE from g^ should dominate, but when 
the numbers of protons and neutrons are equal, N = Z, spin/isospin selection rules cause 
isoscalar ^ interactions, such as go and Ci^ 2 , not to contribute 143 at LO. N = Z nuclei like 


the deuteron are mainly sensitive to isovector couplings, in particular gi, and the EDM from 


9 is suppressed by with respect to the naive expectation in Eq. (52). Using Eq. 

(18) we obtain the order-of-magnitude estimates for dd^t^h in Tab. 6.2, which indicates that 
light-nuclear EDM measurements would offer clues regarding X sources. Expectations about 
the deuteron MQM are described in Ref. 


139 
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In principle, light-nuclear EDMs can be calculated consistently within EFT. Light nuclei 
are sufficiently diluted to be studied in an EFT where pions are integrated out [^, but in this 
case one cannot easily keep track of the chiral-symmetry constraints discussed in Section 
For such nuclei, the pions present in Chiral EFT can be treated in perturbation theory |144| , 
in which case only contact interactions need to be included in the PT sector at LO. For 
the deuteron, the PT 145 and 146 electromagnetic FFs have been calculated quite 
successfully in this framework. This framework was extended to X FFs in Ref. 135 , which 
gives results similar to an earlier calculation based on a zero-range model 147 , and in Ref. 


140 


Similar calculations could be performed for helion and triton. Unfortunately, however, 
pions become non-perturbative at momenta Q ~ TV 148 and thus this approach will fail 
for sufficiently dense nuclei. Much work has been done with nonperturbative pions following 
Refs. 


66 68 , as reviewed for example in Refs. 39,40 . By now, good potentials exist at N^LO 


and N^LO (for example. Ref. 149 ), but for only narrow ranges of UV regulators. Although 


not consistent from an EFT perspective because they lack the necessary counterterms at 
each order and generate amplitudes that are not properly renormalized 41 , these potentials 


produce results in few-body systems that are not very different from phenomenological two- 
plus three-body potentials such as Avl8 150 plus UIX 151 . For a calculation of PT FFs for 


A = 2, 3 systems within this approach, see Ref. 
are discussed in Ref. 


42 


152 . Renormalization issues with currents 


In Tab. 6.2 we give a sample of the most recent evaluations of dd,t,h in terms of the 


nucleon EDMs, Eqs. (38) and (39), the pion-nucleon couplings in Eq. (24), where we bury 


some A contributions, and the other couplings in Eq. (21). For the deuteron, the one-body 
contribution is given by the isoscalar nucleon EDM, 2do = dn + dp^ Isoscalar X potentials, 
generated by go and Ci^ 2 , vanish on the deuteron. go does contribute to dd but only through 
subleading potentials involving PT isospin-breaking couplings, and through two-nucleon X 


currents; these contributions are small 87,136 and omitted here. The gi contribution to dd 
shows little dependence on VpT (see also Ref. 133]), and it is about a factor of 5 smaller than 


the expectation in Eq. (52). A enters only indirectly through gi. 


The situation is strikingly different for helion and triton. The one-body contributions 
are given mostly by, respectively, dn and dpj^ go and gi contribute at about the same level. 
In particular, gi contributes to the isoscalar combination dt + dh, while go to the isovector, 
dt — dh- There is a factor-of-2 disagreement between calculations based on the No-Core Shell 
Model 132 136] and on the Faddeev equation 87,134 . As for the deuteron, the EDM induced 
by the OPE ^ potential is a few times smaller than the expectation in Eq. (52). Changing 


VpT has little effect on the contribution of gi, while go is more affected. The isoscalar couplings 
Cl ,2 give a nonvanishing contribution to dt — dh- These operators are the most sensitive to the 
choice of Vpp and to the details of the nuclear calculation, like the choice of regulator 87,136 , 
and more conclusive results have to wait for a more consistent approach. Still, also this 


^The deviation of the proportionality factor from 1 seen in Table 


6.2 


for the Avl8 and N^LO potentials 

indicating 


stems solely from the deuteron U-state probability Pp. However, Pp is not an observable 153 
that the ~ 10% difference between these potentials and the perturbative-pion result, where Pp enters only at 
N^LO together with other contributions, is within the theoretical error of both approaches. 


^In analogy to the isoscalar trinucleon MDM 153 , one expects model-dependent contributions from the 
trinucleon H-state probability to dt + dh of similar size as those from Pp to dd, see previous footnote. 
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Potential (references) 

dn 

dp 

9l/F-K 

C^Fl 


A/FTrniN 


Perturbative pion 

135 

147 


1 

1 

— -0.23 

— 

— 

— 

dd 

Avis 87, 

131 

,13( 

3 13 

8 

0.91 

0.91 

— -0.19 

— 

— 

— 


NTO 

S'- 

M37 


0.94 

0.94 

— -0.18 

— 

— 

— 


Avis 

132,; 

136, 

138| 


-0.05 

0.90 

0.15 -0.28 

0.01 

-0.02 

n/a 

dt 

AvlS+UlA 

87 

134 


-0.05 

0.90 

0.07 -0.14 

0.002 

-0.005 

0.02 


N^LO 

87 



-0.03 

0.92 

0.11 -0.14 

0.05 

-0.10 

0.02 


Avl8 

132,; 

36, 

138| 


0.88 

-0.05 

-0.15 -0.28 

-0.01 

0.02 

n/a 

dh 

AvlS+UIX 

87 

134 


0.88 

-0.05 

-0.07 -0.14 

-0.002 

0.005 

0.02 


N^LO 

87 



0.90 

-0.03 

-0.11 -0.14 

-0.05 

0.11 

0.02 


Table 3: Dependence of the deuteron, triton and helion EDMs on ^ LECs for various PT 
potentials. Entries are dimensionless in the hrst two columns and in units of efm in the 
remaining columns. “—” indicates very small numbers. 


contribution appears to be smaller than predicted by Eq. (52). A now enters explicitly 


through the three-nucleon^force, Eq. (55), but its contribution is smaller than expected by 
power counting and can be neglected 


Table 6.2| allows us to draw more precise conclusions about the relative size of EDMs, 
which also qualify some of the expectations in Tab. The expectation of relatively large 
light-nuclear EDMs for 9, qCEDM and LR4QOs is tampered by the small effects of the OPE 


^ potential shown in Tab. 6A Still, the even smaller effects from the rest of the P' potential 
leave some opportunities open: 

• 6\ NDA suggests that Ci ^2 and explicit A effects can be neglected. Using the estimates 


for gi in Eq. (27) and the LQCD evaluation of dn^p in Eq. (46), 


dd ~ 0.9{dn + dp) - 0 . 2 ^, 


df dh r\ f\dn dp Qi 


~ 0.9- 


-o-if' 

-*■ 'TT 


(56) 


which satisfy, within large uncertainties, the NDA expectation that the deuteron EDM 
and the helion/triton isoscalar EDM combination receive contributions of similar size 
from the isoscalar nucleon EDM and gi. When dn -|- dp is known, from either experiment 
or a more precise LQCD extraction, a measurement of either dd or dt -|- dh will allow a 


determination of gi. A value of gi much larger than the estimate in Eq. (27) or a dd 


much larger than dn or dp will point to ^ of non-^ origin. The helion/triton isovector 
combination, 

d/f djfi d/fi djn g / ^.\ 

- - h~o.9^ -(57) 

2 2 Ft, 

is expected to be dominated by OPE, but the smallness of g^ and the relative suppression 
of the OPE contribution conspire to enhance the importance of the one-body contribu¬ 
tion. Nonetheless, once dp — dn is measured, Eq. (57) allows the extraction of go, and, 
through Eq. (26), the determination of 6, without further nonperturbative input. 
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qCEDM: As for 6, Ci ^2 and A should be higher order, contributing at N^LO and N^LO 


respectively, so Eqs. (56) and (57) also hold. In this case, however, the couplings are not 
well determined. If we assume that do.i are approximated by their non-analytic pieces 
and that go ^ gi, then dd and dt + dh would be dominated by ^i, with do contributing at 
the 10% level. Due to the logarithmic enhancement in dp—dn and the relative suppression 
of OPE in dt^h) the one-body and OPE contributions to dt — dh should be roughly of the 
same size. Obviously, this is not a hrm conclusion and a better grip on dn,p and ^o,i is 
needed from LQCD. 

LR4QOs: Similar considerations hold, except that the coupling g^ should be small and 
only dd and dt + dh could be enhanced with respect to dn,p. 

gCEDM and PS4QOs: The short-range contributions from (7i 2 are expected to be as 
important as the one-body contribution and OPE, but for the systems considered here 
the effects of the potential are numerically smaller than expected, although highly 
dependent on the short-distance treatment of Vpt- If this result stands, Eqs. (56) and 


(57) hold again and also in this case LQCD input is sorely needed. 


qEDM: the simplest and most predictive situation, in which all EDMs should be well 
approximated by one-body contributions, and both ^ 0,1 can be dropped in Eqs. (56) and 


(57). 


Our discussion has underlined how the EDMs of light nuclei complement the nucleon EDM, 
and can play an important role in singling out the microscopic source of^—for specihc BSM 
examples, see Ref. 154 . Qualitatively, EDMs of nuclei with N = Z are particularly interest¬ 


ing: an observation in these systems of a large EDM, compared to the nucleon EDM, would 
suggest that T does not come from 9, but rather from isospin-breaking sources, like qCEDM 
or LR4QOs. Quantitative conclusions require more precise evaluations of the couplings in Eq. 


(21), a program under way using LQCD. 


7 Outlook 

We reviewed here the main steps that link X in the SM and beyond to observables where 
dramatic experimental progress is expected, nucleon and light-nuclear EDMs. The key idea 
is that EFTs enable us to keep track of symmetries across scales, with minimal assumptions 
about unknown dynamics. As a consequence, the various spin, isospin, and spatial prohles of 
light nuclei probe different aspects of (B)SM and various (admittedly difficult) precision 
measurements would allow precious information on physics at a scale comparable, or perhaps 
beyond, what can be reached in the energy frontier. 

We sketched this link in broad terms, with an emphasis on progress over the last hve- 
ten years. Although the framework is in place, we emphasized a few of the aspects where 
gaps remain. Perhaps the most pressing are i) a. more systematic assessment of the effects of 
flavor-changing operators in the running of the more relevant light-quark operators; ii) LQCD 
calculations of the LECs of the T Lagrangian, either directly or making use of chiral symmetry 
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to relate^matrix elements to T quantities; and in) a fully consistent EFT treatment of the 
PT dynamics in nuclei beyond the deuteron. 

The usefulness of this framework does not hinge solely on experiments on light nuclei. 
Atomic and molecular EDMs (in particular, from diamagnetic atoms) also depend to some 
extent on nuclear^ (especially Schiff moments), albeit for much heavier nuclei than we have 
discussed here. From the nuclear physics perspective, here lies the most important challenge: 
to extend the EFT approach beyond light nuclei. The last decade has witnessed extraordinary 
progress in the development of the so-called ah initio methods, which aim to calculate the 
properties of medium-mass nuclei starting from given internucleon forces. Most of this work, 
just like the trinucleon calculations presented here, treats the input, whether inspired by EFT 
or not, as a phenomenological potential where there is no hierarchy of interactions. With 
such a hybrid approach, the T properties of other nuclei could certainly be calculated using as 
input not only a PT potential inspired by EFT, but also the^ interactions in Eq. (21). This 


would be an alternative to existing calculations (reviewed, for example, in Ref. [^) which are 
typically based on three nonderivative pion-nucleon couplings and uncontrolled (but perhaps 
valid) approximations, such as random-phase or mean-field. We hope this review serves to 
stimulate such alternative calculations. 

On a longer time scale, a fully consistent EFT formulation for larger nuclei would be 
desirable. From an ah initio perspective, nuclear EFT interactions cannot be treated as a 
black box, but instead subleading interactions should be included in perturbation theory to 
avoid extraneous regulator dependence. From a more effective perspective, an in-medium 
EFT appropriate for heavy nuclei should be formulated. The many regularities found among 
nuclear properties offer great opportunities for the effective held theorist. 
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